#include <stdio.h>
#include <stdlib.h>
#include "Operation.h"
#include <math.h>

extern int readFile();
extern FILE* yyin;
/*extern double result[];*/
int main (int argc, const char * argv[]) {
	
	int i;
	double c_phi, c_theta;
	double step_phi =M_PI/32, step_theta = M_PI/32;
	const int size = 2080;/*((2*M_PI)/(M_PI/32))*((M_PI)*(M_PI/32))+1;*/
	double retour[1];
	double paramPhi[1];
	double paramTheta[1];
	double val_rayon[2080];
	double tab_phi[2080];
	double tab_theta[2080];
	int compteur = 0;
		
	for (i=0; i<2080; i++) {
		val_rayon[i]=0;
		tab_phi[i]=0;
		tab_theta[i]=0;
	}
	
	if(argc<2)
		return 1;
	/*calcul pour tout les phi et pour tout les theta la valeur du rayon donnee par l'equation*/
	for (c_theta=-(M_PI/2); c_theta<M_PI/2; c_theta+=step_theta) {
		for (c_phi=0; c_phi<2*M_PI; c_phi+=step_phi) {
			yyin = fopen(argv[1], "r");
			if (yyin == NULL) {
				return 1;
			}
			else {
				/*appel a yyparse*/
				paramPhi[0] = c_phi;
				paramTheta[0] = c_theta;
				readFile(retour, paramPhi, paramTheta);			
				val_rayon[compteur]=retour[0];
				tab_phi[compteur]=c_phi;
				tab_theta[compteur]=c_theta;
				fclose(yyin);
				compteur++;
			}
			
		}
	}
	/*on a nos tableaux de points phi, theta, r*/
	
	
	
    return 0;
}
